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DETERMINATION OF ANGLES OF ATTACK AND SIDESLIP FROM 
RADAR DATA AND A ROLL- STABILIZED PLATFORM 


By John S. Preisser 
Langley Research Center 

SUMMARY 

Equations for angles of attack and sideslip relative to both a rolling and nonrolling 
body axis system are derived for a flight vehicle for which radar and gyroscopic -attitude 
data are available. The method is limited, however, to application where a flat, non- 
rotating earth may be assumed. The gyro considered measures attitude relative to an 
inertial reference in an Euler angle sequence. In particular, a pitch, yaw, and roll 
sequence is used as an example in the derivation. Sample calculations based on flight 
data are presented to illustrate the method. Results obtained with the present gyro 
method are compared with another technique that uses onboard- camera data. 

INTRODUCTION 

Roll- stabilized platforms have been used for many years for vehicle -attitude deter- 
mination on various flight projects. Recently, a miniature attitude -reference system has 
found use in several flight programs at NASA Langley Research Center. (See ref. 1.) 

The system consists of two two-degree-of-freedom gyroscopes which are mounted on a 
common gimbal. The gimbal is roll stabilized; thus, the platform is isolated from the 
rolling motion of the vehicle. An inertial axis system is set up at the time of gyro uncag- 
ing. The orientation of the body axes relative to the inertial axes at any time during the 
flight is determined from the angular displacements measured by the platform after 
uncaging in the sequence pitch, yaw, and roll. 

In order to properly analyze flight-test data for obtaining aerodynamic parameters, 
such as force and moment coefficients and stability derivatives, it is necessary to know 
the orientation of the flight vehicle to the airstream. The purpose of this paper is to 
obtain the equations needed to determine vehicle orientation relative to the airstream 
using an inertial- reference gyro platform in conjunction with radar and wind data. Total 
angle of attack between the longitudinal vehicle axis and the vehicle velocity vector as 
well as the angle between a vector perpendicular to the longitudinal axis and the velocity 
vector has been calculated before from radar data and gyro-platform data using a 
direction-cosine approach. (See ref. 2.) The present method, using an Euler angle 



approach, facilitates the correction of all gy;ro errors, misalinements, and drift prior to 
vehicle lift-off and also provides the components of the total angle of attack — angles of 
attack and sideslip. Angles of attack and sideslip are calculated for both a nonrolling 
and a rolling body axis system. The method, however, is limited to applications where 
a flat, nonrotating earth may be assumed. 

SYMBOLS 

B transformation matrix between earth-fixed axis system and lift-off body axis 

system 

C transformation matrix between earth-fixed axis system and free -stream- 

velocity axis system 

J transformation matrix between inertial axis system and lift-off body axis 

system 

K transformation matrix between inertial axis system and body axis system 

L,L f transformation matrices between free-stream-velocity axis system and body 

axis system 

u,v,w components of free-stream velocity along X-, Y-, and Z-axis, respectively 

u r ,v',w r components of free-stream velocity along X T -, Y f -, and Z r -axis, respectively 

(u f = u) 

Uf,Vf,Wf components 'of free-stream velocity along Xj-, Y^-, and Zf-axis, respectively 

V free-stream velocity (vehicle velocity relative to the surrounding air; see 

appendix) 

X,Y,Z body axis system (a right-handed, orthogonal axis system); X is along the 
longitudinal axis of the vehicle; Y is coincident with the gyro pitch axis 
at the time of gyro uncaging; Z is coincident with the gyro yaw axis and 
in the vertical plane at the time of gyro uncaging 

X r ,Y',Z f nonrolling body axis system; X’ = X; Y T forms a right-handed, orthogonal 
axis system with X r and Z ’ ; the X T Z f -plane is the vertical plane when 

i// f = 0 
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Xj ,Yf ,Zj earth-fixed axis system; Xf is north; Yf is west; Zf is up 

X W ,Y W ,Z W axis system alined with free- stream velocity; X w is along the velocity 
vector; Y w forms a right-handed, orthogonal axis system with X w 
and Z w ; Z w is in the vertical plane with a component up along the 
local vertical 


a 

angle of 

OL ' 

angle of 

P 

sideslip 

P 

sideslip 

t r 
^.o^e.o 

azimuth 


attack, a = tan - *(w/u) 

attack relative to nonrolling body axes, a T - tan _ l(w'/u f ) 
angle, /3 = sin"^(v/V) 

angle relative to nonrolling body axes, y3 r = sin _ l(v’/V) 
and elevation launch angles (launcher settings) 




free- stream flight-path angle in pitch, yp = tan 


-1 


i w f 


/( u f 2 + v £ 2 ) 


1/2 


yy free-stream flight-path angle in yaw, yy = tan-l^-Vf/u^ 

gyro-platform Euler angles measured from inertial reference in the sequence 

O & D 

pitch, yaw, and roll, respectively 

< p,ip,6 Euler angles measured from free-stream-velocity axes to vehicle-body axes 

in the sequence roll, yaw, and pitch, respectively 


Euler angles measured from free-stream-velocity axes to vehicle-body axes 
in the sequence yaw, pitch, and roll, respectively 


Subscripts: 

o at vehicle lift-off 

u at gyro uncaging 
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Matrix notation: 


□ 

O 


square matrix 
column matrix 


A dot over a symbol indicates differentiation with respect to time. 


ANALYSIS 


Angles of attack and sideslip for a vehicle in flight can be derived from a knowledge 
of time histories of range, azimuth, and elevation of the vehicle, the local wind velocity, 
and the vehicle attitude relative to some defined inertial reference. In addition, these 
angles can be derived free from any gyro errors that occur prior to the start of the flight 
test. It will be shown that to accomplish this requires defining five separate axis systems 
and five axis transformations. The axis systems are an earth-fixed axis system, the gyro 
inertial axis system, the flight-vehicle axis system at lift-off, an axis system defined by 
the free-stream-velocity vector, and a body axis system. All axis systems used in this 
report are right handed. 


Transformation 1 

The transformation equations relating the lift-off attitude of the flight vehicle to an 
earth-fixed axis system are 



r ~\ 
X 


k o 




l Z 


°J 


where 



cos r a>0 

sin r a,o 

0 


cos r Gj0 

0 

-sin r Gj0 

B = 

-sin r a?0 

cos r a>0 

0 


0 

1 

0 


0 

0 

1 


sin r 6j0 

0 

cos r e o 
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Multiplying the two matrices yields 



cos r a>0 cos r e>0 

sln r a,o 

- cos r a,o sin r e,o 

B = 

-sin r a,o cos r e,o 

cos r a,o 

sin r a,o sin r e,o 


sin T e o 

0 

cos r e,o 


The axis systems and angles are shown in figure 1. The launch attitude of the vehicle is 
such that the X^Q-plane is perpendicular to the horizontal plane of the earth and the 
Y 0 -axis is perpendicular to the XoZ 0 -plane. The earth-fixed axis system is alined so 
that Xf points north, Yf points west, and Zf points up along the local vertical. The 
launch angles of the flight vehicle at lift-off are then r a0 and r e0 . It is assumed 
that there is no roll between the earth-fixed axis system and the body axis system at 
lift-off. 

x 0 



Figure 1.- Earth-fixed axis system and body axis system at vehicle lift-off. 


Transformation 2 

The radar provides range, azimuth, and elevation data from the radar site to the 
vehicle throughout the flight. These data in conjunction with wind data can be easily con- 
verted to the pitch y p and yaw yy flight-path angles. (See appendix.) It is assumed 
herein that the earth is represented as a flat, nonrotating reference frame. (This 
assumption limits the analysis to short-range, short-duration flights.) The free-stream- 
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velocity axis system is related to the earth -fixed axis system by the following 
transformation: 


*1 


Xw] 

Yf 

> =c< 

Yw 

K 


Z w 
^ J 


where the C matrix will have the same form as the B matrix: 


C = 


cos y y cos y p sin y y 


-sin y y cos y p cos y y 
sin y p 0 


-cos y y sin y p 


sin y y sin y p 


cos y r 


The two axis systems and the angles are shown in figure 2. X w is alined along the free* 
stream -velocity vector. Also, by definition, the X W Z W -plane will be vertical (perpendic-- 
ular to the horizontal plane of the earth) at all times. Since the vehicle velocity is chang- 
ing throughout the flight, y y and y p will be functions of time. 


Zf 



West 


North 


Figure 2.- Earth-fixed axis system and free-stream- 
velocity axis system. 
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Transformation 3 

The body axis system at any time during the flight is related to the uncage position;, 
or inertial axis system, through the general platform transformation: 



where 


K = 


1 

0 


0 


0 

COS 0g 

-sin 0 g 


0 

sin 0g 
cos <£g 


cos Ipg 
-sin tf/ g 
0 


sin ipg 

0 


COS 0g 

0 

-sin fig 

cos l//g 

0 


0 

1 

0 

0 

1 


sin 0g 

0 

COS 0g 


Multiplying the three matrices yields 

-cos t/'g sin 6g 

cos <pg sin i//g sin 0 g + sin 0 g cos 0 g 
-sin <pg sin 0 g sin ipg + cos 4 > g cos 


K = 


cos if/g cos Bg sin i// g 

-cos <pg sin ipg cos 0 g + sin <pg sin Bg cos i// g cos <£ g 
sin <pg sin t^ g cos Bg + cos <£ g sin Bg -sin <f>g cos \pg 


where 0g, and $g are the pitch, yaw, and roll Euler angles, respectively, as 
obtained from the output of the particular gyro platform being considered in this analysis: 
The angles and axis systems are shown in figure 3. 



Figure 3*- Inertial axis system and body axis system. 
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Transformation 4 

For various reasons, a gyro platform might not have zero readings on all three 
axes at the time of vehicle lift-off. For some mission requirements, the gyro might hav$ 
been uncaged while the vehicle was in one position and launched when the vehicle was in , 
another. (See ref. 3.) In addition, since a roll-stabilized platform provides a true iner- 
tial reference in space, the earth’s rotation will cause readings on the gyro while the 
vehicle is still on the earth. The magnitude of these readings will depend on any time 
delay between uncaging and lift-off. Also, residual gyro readings due to gimbal misalinen 
ment and uncage errors may occur. From whatever source, the initial readings can be 
accounted for simply by referencing them back to some effective zero position. If the 
initial gyro readings at lift-off are 0 g>o , and 0 gjO , the lift-off position may be 

considered to have been arrived at by going through the Euler angle sequence 0g )O > ^g,p> 
and (pg o that originated at some effective zero position. The effective zero position is 
the true inertial axis system. The transformation equations can be written as 


r *\ 

Xo 

< Y °> 

Z 0 

v. J 



Z u 


where the J matrix has the same form as the K matrix. (It is assumed that the gyro 
platform is mounted so that it is alined with the vehicle axis.) 



COS q COS q 

sin ^ g)0 

-cos 4, gfQ sin e g>0 

J = 

-cos <£> g>0 sin ^ g> o cos d g>0 + sin <p g)Q sin 0 g>o 

cos ^g,o cos 0g,o 

cos <p gfi sin * gj0 sin 6 gfi + sin <£ g>0 cos 0 g>o 


sin 4>g }0 sin i// g>0 cos 0 go + cos <p gf0 sin 0 g>o 

-sin 0 gjO cos i// g o 

— sin q sin i//g q sin fig q + cos 0g q cos fig q 


Now, by using the four sets of matrix- transformation equations just defined, the 
free-stream-velocity axis system can be directly related to the body axis system. That 
is, 



= KJ' 1 B' 1 



r x^ 

A w 


= KJ -1 B -1 C < 


Y w > 



l Zw J 


where J"* and B~* are the inverse matrices of J and B, respectively. (For 
orthogonal transformations, such as those being considered here, the inverse matrix is 
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identical to the transpose matrix and can be obtained by interchanging the rows and the 
columns of the matrix array.) 

Transformation 5 

A new set of Euler angles can always be defined. If the free- stream- velocity axis 
system is related to the body axis system by a conveniently chosen set of three angles, 
a new transformation equation may be written. For example, 



where L', for the moment, is some new transformation matrix that is a function of the 
newly defined Euler angles. From the previous considerations, 

L’ = KJ" 1 B" 1 C 


Since matrices are equal if and only if the corresponding elements are equal, the elements 
of L f can be expressed in the following manner: 

= Z Kim 

J m 

where k, m, and n will run from 1 to 3 for any i and j. For a given flight applica- 
tion, the J T s and B ? s are constants and need only be evaluated once. The K’s and 
C f s are functions of time. Two convenient sets of Euler angles that can be used to 
define the L f matrix are (a) a yaw, pitch, and roll sequence and (b) a roll, yaw, and 
pitch sequence. 

Transformation (a) .- If the new set of Euler angles were chosen so that the free- 
stream-velocity axes were related to the body axes by a yaw xp\ pitch 0 f , and roll 0’ 
sequence, then the transformation matrix would be as follows: 



L' = 


cos 0 1 cos \p' 

sin 0 f sin 0’ cos \p' - cos 0 r sin \p' 
cos 0’ sin 0 f cos \p' + sin 0 T sin xp' 


cos Q' sin \p' 

sin 0 ’ sin d' sin xf/' + cos c p' cos ip' 
cos 0 f sin 0’ sin ip' - sin <p r cos xp J 


-sin 0’ 
sin 0' cos 0* 
cos 0’ cos 0’ 
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Now by inspection of the L f matrix, ft can be seen that 
^tan-^/L^ 
e' = sm 

cf ) 1 - tan 1 ^2 3^133! 

where ifr' , 6 \ and (/>’ are all functions of time, since the L f matrix is a function of 
time. These angles are defined as Euler angles; however, they can easily be related to 
an angle of attack and sideslip relative to a nonrolling set of axes on the flight vehicle. 

The nonrolling set of axes is shown on the left-hand side of figure 4. These axes 
are an intermediate set obtained during the transformation from the X w ,Yw,Zw set 
to the X,Y,Z set. They are obtained from the free-stream-velocity axes by perform- 
ing the first two rotations only (that is, x ). / and 0 T ). Also shown in the figure are an 
angle of attack c? and a sideslip angle {S' relative to the nonrolling system. Angles 
of attack a' and sideslip /3’ are defined in terms of the vehicle-velocity components 
along the nonrolling axes (see symbol list). However, they are directly related to the 
Euler angles. From the figure, it can be deduced that 

a' = 6 ' and /3 T = -i// 



Figure Free-stream-velocity axis, nonrolling body axis system, 

and body axis system. 
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Finally, the body axes are obtained by rolling an amount, 0 f , from the nonrolling set. 
This axis system is shown in the right-hand side of figure 4. 

Transformation (b) .- By using a transformation different than the previous one, an 
a and a p measured from the body axes can be obtained. The desired angles can be 
obtained by using a roll, yaw, and pitch sequence from the free-stream-velocity axes to 
the body axes. The transformation equations can be written as 



where L now takes the form 

cos xp cos 8 
L = -sin xp 
cos xp sin 6 

By inspection of the L matrix, it can be seen that 
^ = sin _1 (-L 21 ) 

6 = tan'^Lj^Lj,) 

<#> = tan 1 ( l 23/ L 22) 

Now L = L f = KJ - ^B~^C, so that Ljj = Ljj and the formula previously derived for Ljj 
can be used again here. 

Figure 5 presents a sketch of the angles and introduces an angle of attack a and a 
sideslip angle /3. In the sketch, it is assumed that the roll operation has already been 
performed; hence, just yaw xp and pitch 8 are shown. It can be seen that 

a = 8 and p = -\p 

where a is the angle between the longitudinal body axis and the projection of the free- 
, stream velocity V on the XZ -plane, and p is the angle between V and its pro jec- 
tion on the XZ -plane. 


cos (p sin cos 8 + sin <fi sin 8 sin <p sin \p cos 8 - cos <p sin 8 


cos ip cos (p 


sin <p cos \p 


cos (p sin t p sin 6 - sin 0 cos 8 sin <p sin 8 sin t p + cos <p cos 6 


111 



z 



Figure 5*- Free -stream-velocity axis 
and. body axis system. 

It should be noted that the L^j equation is general and therefore can be applied to 
any other gyro system with any Euler angle set. All that is required is that the general 
platform transformation be known for use in the J and K matrices. The equations 
for angles of attack and sideslip will remain the same. 

APPLICATION OF TECHNIQUE 

To illustrate the technique previously discussed, gyro-platform and radar data from 
a previous flight test (ref. 3) were used to derive angle -of-attack and sideslip time his- 
tories. The flight test involved a parachute that was deployed from a payload at a speci- 
fied Mach number and dynamic pressure approximately 60 seconds after vehicle lift-off. 
The data period covered in this example extends from 3 seconds before deployment to 
6 seconds after deployment. The free -stream flight-path angles in pitch and yaw (azimuth 
and elevation) of the payload velocity vector are presented in figure 6. Azimuth and eleva- 
tion launch angles of the vehicle in the lift-off position are also given. Euler angles from 
gyro -platform data are presented in figure 7. These data were obtained from a gyro 
which measured payload attitude in a pitch, yaw, and roll sequence. The payload had little 
pitch or yaw motion prior to parachute deployment. Following deployment, however, the 
payload did experience some near-periodic attitude oscillations. Reference 3 states that 
the parachute did not achieve a steady shape during this data period, which most likely 
contributed to the pitch and yaw oscillations. The payload also had a high roll rate prior 
to deployment, but this rate decreased once the parachute was deployed and the inertia 
properties of the system changed. The initial gyro readings at the time of vehicle lift-off 
are also given in the figure. 
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The results of inputting the data from figures 6 and 7 into the equations previously 
discussed are presented in figures 8 and 9. Figure 8 presents a' and j3 r time histo- 
ries for the nonrolling body axis system as well as a and 0 for the body axis system 
which was rolling. As would be expected, the data all oscillate about the zero level. 

The effect of the rolling motion is readily apparent. The cyclic roll appears directly on 
the a and /3 data. The various roll quantities are shown in figure 9. For this figure, 
a 1-second time period was chosen to pronounce the differences. The quantity <£ was 
not plotted since it would have been nearly coincident with the <$' data for this present 
example involving small angles. 






-3 -2 -1012 3 4 5 6 


Time from deployment, sec 


Figure 8.- Angle-of-attack and sideslip histories. 
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Time from deployment, sec 
Figure 9*- Roll histories. 

Data from a similar parachute flight (ref. 4) were also used to determine pitch and 
yaw time histories. The interval chosen was one which corresponded to vertical descent 
of the parachute system. Pitch and yaw were also determined for this same interval by 
a different method (ref. 5) which made use of film data from a camera mounted near the 
nose of the descending payload. Pitch 6' and yaw \p' from both methods are presented 
in figure 10. In general, there is agreement between the two sets of data. The gyro- 
platform data appear continuous and relatively smooth, whereas the camera data are pre- 
sented by data points (representing discrete film frames) and have some scatter. This 
result is to be expected. Errors in the gyro data would show up mainly as a shift in the 
zero level and would only affect the absolute values of the entire set of data. Relative 
values, over a half- cycle oscillation, for example, would be known to a high degree of 
accuracy. In contrast, each frame of the camera film is read independently; therefore, 
some data scatter, or noise, results. An additional observation to be made from figure 10 
is that the magnitude of the camera data for both positive and negative angles appears to 
be greater than that of the gyro -platform data. This difference implies a calibration 
error in either one or both of the measuring systems. Since the camera was not origi- 
nally intended to be used as a metric instrument and the gyro platform had undergone 
extensive calibration, it is believed that an even closer agreement between the data could 
have been obtained through calibration of the camera system. 

ACCURACY 

The uncertainties that arise in the determination of angles of attack and sideslip 
result from the uncertainties in defining the body axes and the free-stream-velocity axes. 
The orientation of the body axes is determined from the gyro-platform data. The errors 
in platform data are functions of time, accelerations applied to the vehicle, and the vehi- 
cle attitude. The orientation of the free-stream-velocity axes is determined from radar 
data and wind data. Radar errors are functions of range, elevation, and azimuth and the 
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05 


e', deg 


Gyro-platform data (present method! 




Time from deployment, sec 

Figure 10.- Time histories of 9’ and \jr' during parachute descent. 


time rate (jf change of these parameters. Hence, an error analysis must proceed from a 
given flight profile, where the body motions and radar parameters are known. 

An unpublished detailed error analysis for an attitude -reference system like the 
one described in this report has been performed. For the parachute flight tests pre- 
viously described, the results indicate errors in gyro pitch, yaw, and roll to be between 
1° and 2° for the time periods shown in figures 6 to 10. In addition to these gyro errors, 
there is an uncertainty of 0.25 deg/min because the earth's rotation has not been accounted 
for. Radar errors have been estimated from references 3 and 4 to be approximately 0.5° 
in y>y and yp. Hence, the calculated angles of attack and sideslip are estimated to be 
known within ±3° for the examples illustrated in this report. 

CONCLUDING REMARKS 

Equations have been derived to determine angle of attack and sideslip for both a 
nonrolling and rolling axis system for a flight vehicle which is radar tracked and which 
carries onboard a gyroscopic -attitude system. The method is limited, however, to appli- 
cations where a flat, nonrotating earth may be assumed. Specifically, a pitch, yaw, roll 
Euler angle gyro was chosen as the example to be worked out. However, the results can 
easily be applied to other gyro systems if the general platform transformation is known. 
Sample calculations were made, and the results compared well with an independent method 
of attitude determination. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., February 14, 1972. 
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APPENDIX 


CONVERSION OF RADAR DATA TO TRAJECTORY ANGLES 

Typically, a radar determines the position of a vehicle in flight in terms of the slanjt 
range R, elevation angle E, and azimuth angle A from the radar site. The radar site; 
is located some distance (x£ 0 ,y| away from the launcher, where the origin of the eartht 
fixed axis system ^Xf,Yf,Zf^ is located. The radar parameters and the relationship 
between the launcher and the radar site are shown in figure 11. The components of dis- 
tance (xf ?yf »Zf) along the earth-fixed axes between the vehicle and the origin can be 
obtained by transforming the radar spherical coordinates to Cartesian coordinates: 

Xj = R cos E cos A + Xj Q 

y f = -R cos E sin A + yj Q 

C 

Zf = R sin E 


Yf West 


Xf ' 

North 

Figure 11.- Radar parameters and relation of radar site to launcher. 

The components of the vehicle velocity Vr in the earth-fixed axis system are 
found by differentiating with respect to time, 

xj = -RA cos E sin A + R cos E cos A - RE sin E cos A 
yj = -RA. cos E cos A - R cos E sin A + RE sin E sin A 
Zf = R sin E + RE cos E 
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APPENDIX — Concluded 


The components of the free- stream velocity, V, in the earth-fixed axis system can 
be found by combining these last equations with the wind components, 

Uf = Xf - V WjN 
v f = yf - v w,W 

w f = Zf 

where V Wj n is ^h e wind -velocity component from the north and * s the component 

from the west. 

The flight-path angles for the free- stream-velocity vector are now defined in terms 
of these velocity components: 

Y y = tan-^-Vf/uj) 
y p = tan 



The vehicle velocity, wind velocity, and resultant free-stream-velocity vector are shown 
in figure 12 along with the flight-path angles* 



North 

Figure 12.- Vehicle, wind, and free-stream velocities and flight-path angles. 
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